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Abstract 

I find a class of black hole solutions to a (3+1) dimensional theory gravity coupled 
to abelian gauge fields with negative cosmological constant that has been proposed as 
the dual theory to a Lifshitz theory describing critical phenomena in (2+1) dimensions. 
These black holes are all asymptotic to a Lifshitz fixed point geometry and depend on 
a single parameter that determines both their area (or size) and their charge. Most 
of the solutions are obtained numerically, but an exact solution is also obtained for a 
particular value of this parameter. The thermodynamic behaviour of large black holes 
is almost the same regardless of genus, but differs considerably for small black holes. 
Screening behaviour is exhibited in the dual theory for any genus, but the critical 
length at which it sets in is genus-dependent for small black holes. 
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1 Introduction 



The scope of concepts and applications of holography has blossomed in recent years. Real- 
ized in terms of the AdS/CFT correspondence, holography asserts that gravitational dynam- 
ics in an asymptotically Anti de Sitter (AdS) spacetime can be mapped onto a (relativistic) 
conformal field theory (CFT) in one less dimension [1]. This duality conjecture has proven to 
be a useful tool in understanding the behaviour of strongly interacting field theories, whose 
dual description is in terms of weakly coupled gravitational dynamics in a bulk spacetime 
that is asymptotically AdS. 

Over the years holographic duality has been extended beyond high energy physics to 
a much broader class of spacetimes and dual physical systems. Some investigations have 
extended these ideas to asymptotically de Sitter [2] and asymptotically fiat spacetimes [3], 
where holographic renormalization has been shown to be a fruitful tool for understanding 
conserved quantities and gravitational thermodynamics. More recently holography has been 
qualitatively employed in improving our understanding of transport properties of QCD 
quark-gluon plasmas [1] and in the construction of gravity models that are conjectured to 
be dual to various systems in atomic physics and condensed matter physics. For example 
holography has been shown to be applicable to a class of strongly correlated electron and 
atomic systems that exhibit relativistic dispersion relations, and whose dynamics near a 
critical point is well described by a relativistic CFT. The AdS/CMat correspondence is 
being actively employed to study superconductivity |5], the quantum Hall effect [6], and a 
number of other condensed matter systems that can be described by CFTs [7]. 

A key feature of interest is the scaling property 

X ^ Ax (1.1) 

exhibited by fixed points governing the behaviour of various condensed matter systems. 
This behaviour differs from that which arises in the conformal group 

Xt X ^ Ax (1.2) 

and is concretely realized in the (2+1) dimensional theory 

S = [ dt(fx - K(V0)2) (1.3) 



known as the Lifshitz theory. The action (1.3) is invariant under the scaling relation (1.1) 



with z = 2 and has been used to model quantum critical behavior in strongly correlated 
electron systems [HI E] . 



2 



From a holographic perspective the conjectured dual to this system is a gravitational 
theory [10] whose equations of motion yielded solutions with spacetime metrics asymptotic 
to the form 

ds^ = e (^-r'^dt^ + ^ + ^'^x^^ (1.4) 
where the coordinates (t,r, are dimensionless. These obey the scaling relations 

\H A-V X ^ Ax (1.5) 



and the only length scale in the geometry is i. Metrics asymptotic to (1.4) can be generated 
as solutions to the equations of motion that follow from the action 

S = J d^^g (^R-2A- \f^uF^'' - ^H^.rH^"^ - -^e^-""^ B^.F^f,^ (1.6) 

where e^'^"'^ is the Levi-Civita tensor density and F^j, = d[fj_A^] and iJ^j.^ = di^B^r] are 
Abelian gauge fields that are topologically coupled with coupling constant C. The quantity 



A is the cosmological constant. In order to obtain the asymptotic behaviour (1.4) these 
constants are given by jlO] 

A = -^^^^ 2z = {Cir (1.7) 

in terms of the length scale i. 

Recently spherically symmetric solutions to the equations of motion that follow from 



(1.6) were obtained for the z = 2. These solutions included a discrete set of solutions known 
as Lifshitz stars as well as a continuous set of black hole solutions having finite temperature 
[TT]. The purpose of this paper is to extend this work to the case of topological black 
holes, namely those whose event horizons are toroidal or of higher genus as a consequence 
of identifications made in the spacetime [12]. 

These black holes differ considerably from one another depending on the size of the 
event horizon, which is uniquely fixed in terms of their electric charge (or vice- versa). Large 
black holes have thermodynamic behaviour that is essentially the same regardless of the 
genus. However small black holes have markedly different genus-dependent thermodynamic 
behaviour. The entropy of genus- 1 black holes scales linearly with the temperature for all 
values of the black hole radius. Higher genus black holes must have a radius > 1/v^, 



with units defined as in eq. (1.4). As their radius decreases, so does their charge, and 



as Th l/"\/5 these black holes approach their extremal AdS counterparts [13] and have 
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the same asymptotic structure. While most solutions must be obtained numerically, it is 
possible to find an exact black hole solution in the higher-genus case. This solution appears 
to correspond to its zero-mass AdS counterpart, though its asymptotic behaviour is given by 
(1.4). Although all these results are for 2; = 2, generalization to higher z should be possible. 

The solutions found here differ from those with obtained from gravity duals of other 
non-relativistic quantum systems for which topological black hole solutions have also 
recently been found |15j. 



2 Lifshitz Asymptotics and Exact Solutions 



The field equations that follow from the action (1.6) are 

V^F^. = -^-^e^,«;3i^^"^ (2.1) 
C 

'^^HilVT = r) r^ ^llualiF"''^ (2.2) 

G^u - ■^Qf.u = 2 (^^f^rF^ ~ ^9iiuF'^^ + - (^H^^rH/^ - -g^uH"^ (2.3) 

In order to solve these equations an ansatz is needed that preserves the basic symmetries 
under consideration. For the metric I shall take 



ds^ = f I -r^f\r)de + + rHnl ] (2.4) 



where 

de"^ + sm^ edcf^ k = +l 

dnl = { de^ + e^d(j)^ k = o (2.5) 

d9^ + sinh^ edd)^ k = -1 



is the metric for spatial sections at fixed (t,r) corresponding to the genus g = (spherical), 
g = 1 (fiat/toroidal), and g > 2 (hyperbolic) cases. Compact spatial sections at fixed r are 
possible by making appropriate identifications in the x coordinates [16], [17] . 
The gauge field strengths are 

{sin6' k = +l 
6 k = (2.6) 
sinh6' k = —1 

with all other components either vanishing or being given by antisymmetrization. 
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The field strength is that of an electric field directed radially outward, and described 
by the function h{r) in an orthonormal basis. This field sources the 3-form field strength 
H^^T and vice-versa. This latter field is therefore electrically charged and corresponds to a 
charged fluid whose density is given by j{r) in an orthonormal basis [TT] . 



The field equations (2.1 
df 



2.3) reduce to the system 



dr 
dg 
dr 

^dj_ 
dr 
dh 
dr 



Im + lf{r)g{rr 



h{r) 



1 



5+ -^-j{rf - h{ry 



2^ 
2h{r) 



k 



2g{r)h{r) + -j{r) - -j{r)g{r) [5+ — +j{r) - h{r)' 
2g{r)j{r] 



(2.7) 
(2.8) 
(2.9) 
(2.10) 



which is a system of ODEs that can be solved by standard numerical methods. 
An exact solution to the above equations is 



fir) 



yielding the metric 



J[r) 



1 + 



2^ 



1 



h{r) 



(2.11) 



ds' 



r 



k 

£2 + 2 



dt"^ + 



dr'^ 



+ r'dQl 



(2.12) 



Kt^ ~^ 2) 

where I have rescaled r —>■ r/i and t —>■ t/i. This solution is valid only for z = 2 and does 
not straightforwardly generalize to other values of z > 1. 

The spacetime described by (2.12) is singular at r = for all values of k. For k = 
it becomes the metric (1.4): all its curvature scalars are finite but it is not geodesically 
complete . This kind of metric can be regarded as physically reasonable if there exists a 
regular black hole solution that approaches it in some extremal limit [T0l[T8]. For k ^ the 
Kretschmann scalar diverges are r = 0. If = 1 this is a naked singularity, but if /c = — 1 
it is cloaked by an event horizon. Hence the metric 



di' 



dr + 



dr^ 



(9-1) 



(2.13) 



is an exact black hole solution to the field equations provided h{r) = 1 and j{r) = 
The event horizon is located at r = i/\/2. This metric would appear to be the analogue of 
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a zero-mass topological AdS black hole [13], though it should be noted that defining mass 
and other conserved charges in this theory remains an open question at this stage. 

The asymptotic behaviour of the system (2.7)-(2.10) is the same for all values of k, and 
so its analysis is identical to the given for the k = 1 case [H]. Equations (2.8)-(2.10) form 
a closed system for the set {g{r), j{r), h{r)} and can be considered separately; once they 
are solved then eq. (2.7) can be solved for /(r). Linearization of these equations indicates 
that there is a zero mode at large r. However this mode must have zero amplitude if the 
system is to approach the Lifshitz metric (1.4). This can only happen if the initial values 
of the functions are appropriately adjusted (fine-tuned) to ensure that each function in the 
set {g{r),j{r),h{r)} approaches unity. 



3 Series and Numerical Black Hole Solutions 



The system (2.7)-(2.10) can be solved both for large r and near the event horizon using 



series expansions. For large r the solution is 

k k^ / k^ khL\ 1 nSk'^hL hk^ 15/it 



4r2 32r4 \^128 12 y \^ 384 2048 64 J 



4r2 \ 2 32 y r4 \^128 8 / ^ 32 2048 32 y r 



4r2 \ 2 32 J \128 8 J \ 32 2048 32 J 

r4 2r6 V 32 ^ 16 ^ ^ 

and is governed by one constant h^, that parametrizes how rapidly the electric field falls off 
at large r. For k = the falloff is very rapid, particularly for /(r), whose subleading term 
is proportional to 1/r^. 

It is also possible to solve the equations near the event horizon. Under the assumption 
that the black hole is not extremal, the metric functions Qu and grr must respectively have a 
simple zero and a simple pole at the horizon r = tq. The electric field /i(ro) = is assumed 
to be finite at r = tq. Unlike the situation for the (AdS)-Reissner-Nordstrom metric, these 
values are not independent because of the presence of the charged fiuid (the 3-form field 
strength). For a static configuration the gravitational pull of the black hole must balance 
the self-repulsion of the fiuid, whose value must also vanish at the horizon if the electric 
field is finite there. 
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Writing r = tq + x, the resuhant near-horizon series solution is 



where 



/i 

^0 

jo 
jl 



f = /ov^ (1 + M + /2X2 + ■ ■ ■) 

g = (l + gix + g2X^ ^ ) 

y/X ^ 

j = joVx (l + jlX + j2X^ H ) 

h = ho {l + hix + h2x'^ -\ ) (3.2) 

5k^ - 524hl + 45A:rg - Ukr'^hl + 6r^h^ + lOOr^ 
2{ro{k + 5r2 - rg/ig)^ 

3/2 



3A;2 - 2Ar^hl + 2bkrl - 7kr^hl + Ar^h^ + 50rf 







2ro{k + 5rQ — r^h^y 



ho{e + llfcrg - krlhl - Sr^hl + 30r^ 



with the remaining coefficients determined in terms of ho and tq. 



From (3.3) it clear that the constraint 

(5 - hl)rl + k>0 (3.4) 

must be respected so that the solutions are real. For k = 0,1 this imposes an upper bound 
on \ho\, whereas for k = —1 it imposes a lower bound of ro > 1/V5 on the size of the black 
hole. 

Further progress can only be made numerically. Since the equations are all ffist order 
ODEs, the shooting method can be used to solve them. Choosing some value of tq for the 
black hole, one can then choose a value for ho and then find the initial values of the functions 
g{r), h{r), and j{r) from the series solutions for some value of r > tq. This becomes the 



initial data for the system (2.8)-(2.10), which can then be numerically solved. The value of 
h{r) can be computed for r = R » tq; should \h{R) — 1\ > e, where e << 1 is some level 
of tolerance, then the process is repeated, with value of ho adjusted, until \h{R) — 1| < e. 
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Figure 1: ho vs. tq on a log-log scale for spherical {k = 1, solid), toroidal {k = 0, dot) and 
higher-genus {k = —1, dot-dash) black holes. The limit given in eq. (3.4) for = 1 is the 



dashed line. The intersection point of the k = —1 curve with the tq axis corresponds to the 



exact solution (2.13). 



This then gives a black hole solution of radius tq with asymptotic behaviour given by the 



metric (1.4) 



In practice it proves best to begin with large black holes, for which the value of Hq remains 
constant. Numerically I find that setting ho = 1.374 yields valid solutions for tq > 10 for 
all values of k. For k = this value of ho yields black hole solutions for all values of tq that 
were computationally viable. For tq < 10 it is necessary to systematically adjust the value 
of ho upward for k = 1 and downward for k = —1. 

Figure [T] illustrates the behaviour of ho as a function of tq on a log-log plot. For k = 1, 



the value of ho asymptotes to the upper bound given in (3.4). The extremal limit appears 
to occur when the bound is saturated as tq ^ 0; it does not appear to be possible to find 
series solutions for tq > in the extremal case. For k = —1 the value of ho approaches 
zero as tq — > 1/ V5. In this case the limit is evidently that of an extremal topological AdS 

I/V2 and ho = h{r) = 1, 



black hole with genus g > 2. For the exact solution (2.13), ro 
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Figure 2: A plot of the metric and gauge functions for ro = 20 for all three values of k. The 
three curves overlap within the plotting resolution. 

corresponding to the intersection point of the A; = — 1 curve with the ro in figure [T] The 
exact solution provides a useful check on the numerical solution, and it can be verified that 
the solutions for the metric and gauge functions are identical to within limits of tolerance 
(taken to be 10~^). 

The metric functions and gauge field strengths are essentially indistinguishable amongst 
the different values of k for tq > 10, as an example illustrates in figure |2} Small differences 
begin to appear for intermediate values of tq, as shown in figure [3} and for small Tq the 
distinctions are quite significant, as depicted in figures |4] and [5] 



4 Black Hole Thermodynamics 

The temperature of Lifshitz black holes is easily evaluated using standard Wick-rotation 
methods, yielding the result 

4vrc/o 47r 

where /o is dependent upon tq so that the metric has the asymptotic behaviour given in 



equation (1.4). 
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Figure 3: A plot of the metric and gauge functions for tq = 1.6 for A; = 1 (solid), A; = 
(dot), and k = —\ (dot-dash). 
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Figure 4: A plot of the metric functions for tq — 0.4 iox k — 1 (solid), k — Q (dot), and 
k — —1 (dot-dash). 
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Figure 5: A plot of the gauge functions for tq = 0.4 for A; = 1 (solid), A; = (dot), and 
k — —1 (dot-dash). 




Figure 6: Dependence of the black hole temperature as a function of Tq for k = 1 (solid), 
A; = (dot), and k = —1 (dot-dash). 
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Figure 7: Dependence of the black hole entropy as a function of temperature for k = 1 
(solid), k = (dot), and k = —1 (dot-dash). 

The temperature for any genus is therefore dependent only on the size tq of the black 
hole, and its behaviour is illustrated in figure [6] For large tq, spherical black holes {k = 1) 
are hotter than planar /toroidal black holes {k = 0), which are in turn hotter than the 
k = —1 higher genus class. The planar /toroidal holes have a temperature dependence that 
increases quadratically with the black hole radius; numerically I find that 

Tfc=o = 10-°-^Vo' = .276rl (4.2) 

a behaviour that also accurately describes the k = ±1 cases for tq > 0.6. This behaviour is 
consistent with the near-horizon expansion that indicates (70 ~ V/o ~ fo^' large tq. 

For small tq the behaviour is strikingly different for each case. Higher genus black holes 
approach zero temperature at an exponentially rapid rate as tq l/V^- For ro < 0.1 the 
planar black holes become hotter than their k = counterparts. The temperature of the 
k = 1 black holes falls off very rapidly as tq decreases, though less so than the k = —1 case. 

The entropy of the black holes is given by | of the area, so S = vrrg for all cases (assuming 



appropriate identifications and volume normalizations for k = 0, —1). Hence eq. (4.2) gives 



S = 11.4T (4.3) 
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a relationship that is vahd for all T for the k = case and that holds for large T for /c = ±1 
[11] . The behaviour of the entropy as a function of temperature is given in figure [7j 



For the exact solution (|2.13|), the temperature and entropy are 

S 



T 



1 



2 



(4.4) 



which can be computed directly from the exact solution or from the series solution (3.3) for 
which /o = l/5'o = 2^^^ and rg = 1/ \/2. These values are commensurate with those in figure 
ItI though they are considerably beyond the range for which eq. (4.2) applies. 



5 Wilson Loops and the Boundary Dual Theory 

Since in (2+1) dimensions one can write = V x E, where Ej = ejkV^<f), the boundary 



theory (1.3) can be regarded as a gauge theory in (2+1) dimensions (albeit one with an 



unusual action) with a dimensionless coupling constant P, [H]. As in the k = spherical 
case, one can introduce Wilson loops by joining charged particles on the boundary that 
are connected together in the bulk via a string. The Euclidean action of this string for a 
rectangular Wilson loop is the same for all values of k and is given by [TT], Hn] 



5 = ^ y dadT^det[gABd,X^d,XB] 



f de\ Pr^'+^ + Pg^r^'-^ ( — 
27ia' J V ^ \d9 



taking a = 6 and ir = t in the static gauge, with Euclidean time interval A. 
Extremizing the action yields a constant of the motion 



from which can be computed the boundary length 

dr g 



z+l 
m 



L 



de 



22+2 



- 1 



and the regularized potential energy between the two particles 

/ 

S f 

V 



s_ 



27ra' 



„2-l 



2 / dr 



fg 



V 



1 



if) 



= -2 / drr'-'fg 

22+2 Jro 



(5.1) 



(5.2) 



(5.3) 



(5.4) 
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Figure 8: Dependence of the boundary length L between two particles as a function of the 
string midpoint for = 1 (solid), = (dot), and k = —1 (dot-dash), where tq = 0.5. 
These curves rapidly become indistinguishable as tq increases. 



where fm = fijm) and > tq is the location of the midpoint of the string. 

For small the separation L between the particles grows, reaching a maximum and then 
decreasing as gets larger. The energy between the particles will is negative for sufficiently 
small and will vanish at some particular value of (or L). Beyond this point the energy 
for a single string joining the pair is positive. The energy of the configuration is therefore 
minimized (to zero) by two (non-interacting) strings stretching from each particle down to 
the horizon, screening the gauge interaction between the particles. The point at which this 
takes place will depend on ro, and hence the temperature. 

For large black holes the metric functions / and g are nearly indistinguishable for all 
values of fc, and so the analysis of screening behaviour for the spherical = case [TT] holds 
for the other values of k as well. For small black holes the value of at which the potential 
vanishes is numerically almost the same, as in figure [8j However the length between the 
particles differs considerably as a function of (see figure [9]), and so the critical value 
L = Lc at which screening occurs will differ considerably for different values of k. 
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Figure 9: Dependence of the potential V between two particles as a function of the string 
midpoint for A; = 1 (solid), k — ^ (dot), and k — —\ (dot-dash), where ro = 0.5. The 
right-hand diagram is a close-up of the left. The critical values of for each k are within 
0.04% of each other. The intersection points (and the behaviour of V) rapidly become 
indistinguishable as Tq increases. 
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For the exact solution (2.13) it is straightforward to carry out the integration in equations 



(5.3) and (5.4), though the results cannot be obtained as an explicit function of but must 
instead be obtained numerically. It is straightforward to check that the integrands (and 
hence the integrals) in equations (5.3) and (5.4) are identical within limits of tolerance for 
the exact solution (2.13) and its numerical counterpart with tq = 1/v^- The behaviour is 
not too different for this case as compared to its k = 0,1 counterparts, and so the results 
are not plotted here. 



6 Conclusions 

In this paper I have expanded the class of black holes that provide a dual description of a fi- 
nite temperature a Lifshitz system to include black holes of any topology. The gravitational 
theory can at least be regarded as a phenomenological description of the 2+1 dimensional 
physics described by the Lifshitz theory. Whether or not this duality can be fully incor- 
porated into string theory, thereby extending the AdS/CFT correspondence, is an open 
question. 

Thermodynamic properties of these black holes are quite similar for large black holes, 
but differ considerably for small black holes. The genus and 1 cases approach extremality 
as the black hole size approaches zero. However the higher-genus black holes approach 
extremality at ro = 1/ v^, at which point the gauge fields vanish and a negative mass 
topological AdS black hole is attained. 

The screening behaviour of the dual theory is essentially the same for any genus for large 
black holes. For small black holes the onset of screening is attained at the nearly same value 
of the midpoint of the string joining two charged particles on the boundary regardless 
of the genus; however this yields very different values of the critical length L between the 
particles on the boundary due to a sensitive dependence of this quantity on Vm- 

A number of interesting questions remain, including higher-dimensional generalizations, 
a more complete study of the z ^ 2 cases, developing holographic renormalization for this 
class of theories, and understanding better the relationship with non-relativistic, non-abelian 
gauge theories having quantum critical behavior at 2; = 2 [20] . 
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